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Mh Most excited hadrons have multiparticle strong decay modes, which 

fi , can often be described as resulting from intermediate states containing one 

^ or two resonances. In a theoretical approach, such a description in terms 

,S^ of quasi-two-particle initial and final states leads to unitarity violations, 

because of the complex masses of the involved resonances. In the present 

,__! paper, an empirical algebraic procedure is presented to restore unitarity 

^ of the iS-matrix while preserving its symmetry. Preliminary results are 

(^ presented in a first application to S-wave tttt scattering, in the framework 

On of the Resonance-Spectrum Expansion. 
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1. Introduction 



. . Inspection of the Particle Data Group (PDG) tables [1] reveals that most 

^^ excited mesons and baryons have strong decay modes involving more than 

S^ two lighter hadrons. Moreover, one also verifies that many of these decays 

can be considered "cascades" of decays from intermediate states involving 
one or more resonances. A few mesonic examples of such processes are [1] 

1. /o(1370) ^ p/9,2(^7r)s-wave, ... ^ 47r ; 

2. i^2(1770) -^ K^ {U30) TT, Kf 2{1270), ... -^ Ktttt ; 

3. (/)(2170) -^ 0/o(98O) -^ (pTTTT, KKfo{980) -^ KKtttt . 
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In the first example, 2(7r7r)5_„ave stands for aa, wliere a is the very broad 
/o(600) [U scalar resonance. In all three cases, the final state contains 3 or 
even 4 mesons that are stable with respect to strong interactions. In order 
to describe such processes in a mathematically rigorous way, respecting uni- 
tarity, one would have to resort to (relativistic) Faddeev [2j or Yakubovsky 
[3] equations, respectively, in the scattering regime. As a matter of fact, 
such a Faddeev approach was applied [1] to the (/>(2170), as a three-body 
resonance in (pKK and (pinr, though not including the KKfQ(980) compo- 
nent and also not the 4-body KKtttt channel. Note that solving relativistic 
3- or 4-body equations with many channels is already a huge task, but it 
becomes absolutely impracticable if one wants to go beyond an effective 
description as in Ref. [4j, by taking into account the quark substructure of 
the decaying resonance. 

In experimental data analyses, one often falls back upon dispersive or 
purely phenomenological parametrisations, as e.g. in Refs. O [6] for the 
description of the 47r channel in 5- wave tttt scattering above 1 GeV and the 
/o(1370) resonance. 

An alternative, theoretical approach is to interpret the intermediate 
state in a cascade decay as a (quasi-)final state, containing one or two res- 
onances as outgoing particles. The problem then arises how to deal with 
the, in principle, complex mass(es) of the resonance(s), while preserving 
two-body unitarity. One way is by discretising the resonance real-mass dis- 
tribution(s), and accordingly including a large number of channels having 
the corresponding threshold energies, with relative coupling strengths given 
by the mass-distribution function. Something in this spirit has been carried 
out in a coupled-channel description of S'-wave meson-meson scattering in a 
chiral unitary model [7]. A drawback of such a method is the proliferation 
of channels, besides the problem of properly dealing with resonances not far 
above threshold. 

In the present study, we shall try out a novel, intuitive approach, by 
simply substituting complex values for the masses of the final-state reso- 
nances, according to their listed [Ij real (Breit-Wigner) masses and widths. 
Although this step inexorably destroys unitarity, we can afterwards restore 
it by a suitable redefinition of the 5-matrix. This quite general procedure 
is then applied to a concrete test case, in the framework of the Resonance- 
Spectrum Expansion. 

2. Resonance-Spectrum Expansion 

The Resonance-Spectrum Expansion (RSE) is a model for the scattering 
of two mesons in non-exotic channels, via an infinite set of intermediate s- 
channel qq states, i.e., a kind of Regge propagators [8j. The confinement 
spectrum for these bare qq states can, in principle, be chosen freely, but in 
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all successful phenomenological applications so far we have used a harmonic- 
oscillator (HO) spectrum with flavour-indepedent frequency (see Ref. [8] for 
several references) . Because of the separability of the effective meson- meson 
interaction, the RSE model can be solved in closed form. The relevant Born 
and one-loop diagrams are depicted in Fig. [TJ For N meson-meson channels 



W< 



vnv : 




Fig. 1. Born and one-loop term of the RSE effective meson- meson interaction, 
and several qq channels, the effective potential has the form 

a=l n=0 E — En 

^ n,,{E)Jl{p,ro)ji^{p',ro), (1) 

where the RSE propagator lZij{E) contains an infinite tower of s-channel 
bare qq states, corresponding to the spectrum of an, in principle, arbitrary 
confining potential. Here, En' is the discrete energy of the n-th recur- 
rence in qq channel a, gl' (n) [9] is the corresponding coupling to the i-th 
meson-meson channel, j^Xpi) is the Lj-th order spherical Bessel function. 
Pi is the relativistically defined off-shell relative momentum in meson-meson 
channel i, tq is a distance parameter, and A is an overall coupling constant. 
Note that the spherical Bessel function results from assuming ^Pq qq pair 
creation/annihilation only to take place at a certain distance tq [8]. The 
closed-form off-energy-shell 7~-matrix then reads 



-2AVo ,/Jmp^j JlM^o) E '^irn{E) {[1 - ^^]"'}^. Ji^ip'jro) , 



T^-^^>{p„j/^,E)= (2) 



m=l 



where 

niE) = -2aVodiag(i£„(A:„ro)/ig"(A;„ro)) , (3) 

with h]^ {knTo) the spherical Hankel function of the first kind, and kn 
the on-shell relative momentum in meson-meson channel n. Finally, the 
corresponding unitary and symmetric (on-shell) 5-matrix is given by 

sl'^-'-'Xh^kj-E) = 5ij + 2iT^^-'''\ki,kf,E) . (4) 

Note that the 5-matrix is only unitary for real /cj, and so for real meson 
masses in channel i. 
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3. Redefining the «5-niatrix 

As mentioned above, if we take a complex mass in any of the meson- 
meson channels, the 5-matrix ceases to be unitary, but stays symmetric. 
This requires a redefinition of the 5-matrix. Now, note that 

S^S = A (5) 

is not unity anymore, but it is still a Hermitian matrix, by definition, and 
so with real eigenvalues, which moreover are all positive in this special case. 
Thus, A can be diagonalised by a unitary matrix U: 

Ad = US^SU^ . (6) 

Define now 



S = SU^A-^^^^U , (7) 

3cause 
forward to show that 



— 1/2 

where ^^ is real, because of the positive eigenvalues. Then, it is straight- 



S^S = SS^ = 1. (8) 

It is not so easy to prove that S is also symmetric, but this has been checked 
numerically with a precision of better than one part in a trillion. 

4. Preliminary results 

Let us now use the redefined 5-matrix of Eq. ([8]) in a purely comparative 
calculation. Starting point is our [10] application of the RSE to 5-wave tttt 
scattering, with pseudoscalar-pseudoscalar, vector- vector, and scalar-scalar 
channels included. Then, we substitute complex values for the masses of 
the broadest resonances in some of these channels, namely for the /o(600) 
(alias a), i^g(800) (alias k), and p meson. Note that the precise values we 
take for the complex masses of the very broad a and k resonances are not 
so important here, but in future model fits to the data one should use the 
best available, "world- average" values for the corresponding pole positions. 
In Fig. [2| we plot the curves for the two cases. Quite significant differences 
become evident, especially above 1 GeV. Nevertheless, the curve obtained 
with S in Eq. ([7|) and complex masses looks reasonable, though no refit has 
been done. 

5. Summary 

In the present study, we have carried out an empirical procedure in order 
to restore unitarity of a coupled-channel 5-matrix with complex masses in 
the asymptotic states. Although no rigorous mathematical justification is 
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Fig. 2. RSE calculation of S'-wave tttt phase shifts. Red (upper) curve: from 
standard S'-matrix and real a, k, and p masses; blue (lower) curve: from redefined 
iS-matrix in Eq. ([7]) and complex masses. 



given here for the purely algebraic transformation of the original 5-matrix, 
preservation of the mandatory symmetry of S gives us confidence that the 
method makes sense. A direct comparison of S'-wave vrvr phase shifts, first 
calculated with the original S and using real masses for the a, p, and k 
resonances in some of the coupled channels, and then with the S of Eq. (JTl) 
and complex values for the latter masses, shows significant yet reasonable 
changes. Nevertheless, our procedure will have to be tested in concrete fits 
to the data so as to find out whether it is really a promising new approach 
to multiparticle hadronic decays. 
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